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Abstract. We present a new code for calculating the Fe Ka line profiles from 
relativistic accretion disks with finite thickness around a rotating black hole. 
The thin Keplerian accretion disk must become thicker and sub-Keplerian with 
increasing accretion rates. We here embark on, for the first time, a fully rel- 
ativistic computation which is aimed at gaining an insight into the effects of 
geometrical thickness and the sub-Keplerian orbital velocity on the line profiles. 
This code is also well-suited to produce accretion disk images. 



Introduction 



The fluorescent Ka iron emission line has been observed in several active galac- 
tic nuclei(AGN). This line is thought to be produced by iron fluorescence of a 
relatively cold X-ray illuminated material in the innermost parts of the accretion 
disk. Thus the line profile encodes the nature of the structure, geometry, and 
dynamics of the accretion flow, as well as the geometry of the space. 

Calculations of the line profiles have been performed by several authors. 
Two basic approaches have b een described: 1) using a kind of transfer functio n 
to map a disk into sky pla ne (Cunningham 19751; [Laorlll99ll: ISpeith et al.| [l995); 
2) a ray tracing approach ( Cadez et al.lll998l : lBeckwith fc Dondl2004l ) . But all of 
these techniques are restricted to the geometrically thin, Keplerian disk located 
in the equatorial plane. To our knowledge, the disk must bec ome geometri- 
cally thick and be sub-Keplerain wit h increasing accretion rates (jWang &: Zhou 
19991 : IShadmehri Khaienabi 20051 ) . Motivated by the above considerations a 
relativistic thick accretion disk model is presented. 



2. Assumptions and Results 

For simplicity, we adopt a conical surface for the disk geometry. The thickness 
of the disk can be described by the declination angle 6(0 < 5 < vr/4). When 
6 = 0, the disk reduced to a thin disk. For sub-Keplerian velocity, w e adopt the 
modification of firstly introduced by iRuszkowski Sz Fabian ( 20001 ) 
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It is easy to verify that Q < £Ik, is the Keplerian angular velocity. The 
parameters of this model include: the radii of the emitting disk zone r- m , r ou t; the 
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spin of the black hole a; the inclination angle of the disk($ ) and the angle(i? e ) 
between the rotation axis of the system and the radial direction of the disk 
surface which satisfies 5 + i? e = vr/2; the radial emissivity index p; the angular 
velocity index n, respectively. 

Following the idea presented bv lSpeith et al. ( 19951 ). we extend their model 
to the finite thick disk, and adopt elliptic integrals to improve the performance 
of the code. Here, we show, as an example, a number of line profiles obtained by 
our code(see FigfT]), the emissivity law in all cases is taken the forms e(r e ) cx r~ 3 
and f(ne) oc (1 + 2.06/x e )- 




Figure 1. Examples of line profiles showing the impact of different model 
parameters to the profile. The emission line region is from r ms to r ou t = 20r g , 
where r g is the gravitational radius. Top left: line profiles as a function of the 
angular velocity for n = 1,2, 3, 4. Top right: line profiles as a function of the 
disk thickness for tf c = 50°, 60°, 70°, 80°, 90°. Bottom left: Comparison of the 
relativistic line profiles generated by our model with different observer inclina- 
tions i? = 15°, 30°, 45°, 60°, 75°. Bottom right: Comparison of the relativistic 
line profiles generated by our model with different spins a = 0, 0.5, 0.998 and 
n — 3. The flux in all cases is given using the same arbitrary units, and all 
our results are unsmoothcd. 
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